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ABSTRACT 


A  compressor  designed  for  LDA  measurement  of  unsteady  rotor 
flow  has  been  constructed.  A  blower  in  series  uncouples  throughflow 
and  compressor  speed  providing  wide  variability  of  load.  LDA  data 
acquisition  and  processing  systems  have  been  designed  and  partially 
completed.  Candidate  data  processing  algorithms  have  been  studied. 
Propagating  stall  has  been  detected  experimentally  but  no  LDA  data 
have  been  obtained  defining  unsteady  velocities. 

Modelling  of  a  cascade  with  doubly  periodic  flow  has  been 
initiated.  Solutions  with  nonuniform  loading  have  been  demonstrated; 
future  addition  of  the  vortical  flow  downstream  should  lead  to 
an  instability  criterion  and  the  subsequent  development  of  time 
dependence. 

Unsatisfactory  progress  has  prompted  project  reorientation 
with  a  CFD  cascade  model  which  creates  and  convects  vortex  "blobs" 
according  to  the  Random  Vortex  Method.  This  model  may  be  merged 
with  the  previous  analytical  model  of  nonuni formly  loaded  cascades. 

Compressor/LDA  development  will  continue  with  emphasis  shifting 
to  a  windtunnel  experiment  with  one  stalling  airfoil.  Measured  flow 
periodicity  will  validate  the  CFD  model  and  allow  application  to 
cascades.  This  objective  with  3D  application  will  allow  predicting 
the  stall  zone  number  in  compressors,  the  missing  datum  in  assessing 
blade  resonance  with  stall  frequency.  The  tangential  wavelength  is 
also  key  to  predicting  stall  flutter. 
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INTRODUCTION  AND  BACKGROUND 


The  axial-flow  compressor  and  fan  components  of  aircraft 
gas  turbine  engines  are  typically  designed  for  minimum  weight 
and  high  efficiency.  In  the  aircraft  application  frequent 
operations  at  off-design  conditions  are  required,  resulting  in 
the  occurence  of  potentially  dangerous  stall,  or  near  stall,  in 
these  compression  components.  The  danger  is  related  to  the 
possibility  of  vibration,  fatigue  and  eventual  failure  of  the 
axial-flow  blading.  These  vibrations  may  be  self-excited  in 
which  case  they  are  termed  "stalling  flutter",  or  else  they  may 
result  from  a  resonance  with  the  frequency  of  stall  encounter 
in  the  "stall  propagation"  mode.  Hence  the  detailed  mechanism 
of  blade  stall  in  a  cascade  of  airfoils  needs  further  definition. 
Present  understanding  is  deficient  in  many  respects. 

Deficiencies  with  respect  to  stall  propagation  relate  mainly 
to  the  frequency  of  stall  encounter  by  a  blade.  Although  the 
propagation  rate  of  a  stall  patch  along  the  cascade  may  be  pre¬ 
dicted  quite  accurately,  the  number  of  such  patches  in  the  flow 
annulus  (i.e.  wavelength  along  the  cascade)  cannot  presently  be 
predicted.  Hence  the  product  of  rate  and  wavelength  which  de¬ 
termines  the  forcing  frequency  is  not  accessible  to  the  design 
process . 

The  deficiency  with  respect  to  stalling  flutter  is  even 
more  severe;  neither  the  precise  onset  of  the  phenomenon  nor  its 
severity  can  be  predicted  from  first  principles.  It  is  known  to 


be  related  to  operation  at  high  incidence,  in  the  proximity  of 
stall,  and  periodic  stalling  and  unstalling  is  thought  to  take 
place  during  each  cycle  of  oscillation.  Recently  compressibility 
phenomena  have  been  implicated  in  the  stalling  process  in  cer¬ 
tain  instances. 

As  an  overall  statement  it  seems  clear  that  a  better  de¬ 
finition  is  needed  of  how  a  cascaded  airfoil  stalls  dynamically. 
Consequently,  the  research  program  reported  here  has  had  as  its 
primary  objective  the  detailed  measurement  of  the  velocity  field 
in  a  compressor  rotor  passage.  It  is  clear  that  this  is  an  un¬ 
steady  velocity  under  the  condition  of  interest.  Furthermore, 
the  highest  positive  incidence  is  usually  obtained  on  the  first 
rotor  row  of  blades,  hence  stalling  flutter  and  instigation  of 
propagation  stall  may  be  studied,  initially,  in  a  single  stage 
compressor.  Multistaging,  moving  shockwaves,  and  other  effects 
can  be  left  to  future  consideration. 

With  this  background  the  basic  thrust  of  the  research  was 
intended  to  be  experimental.  A  research  compressor,  specifically 
conceived  for  this  purpose,  was  designed  and  built.  A  laser  dop- 
pler  anemometer  was  converted  to  the  service  of  these  turbo¬ 
machinery  measurements,  and  a  data  acquisition,  storage  and 
processing  system  was  procured.  These  efforts  are  detailed  in 
the  following  sections. 

Finally,  in  anticipation  of  obtaining  a  complete  velocity 
mapping  of  a  rotor  passage  in  three  dimensions  and  time,  a  pre¬ 
liminary  modelling  of  the  flow  was  begun  in  order  to  help  to 
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guide  later  experiments  and  assess  the  difficulty  of  the 
desired  analytical  description. 

Many  of  the  preliminary  phases  of  the  research  were 
reported  in  [13#  the  final  technical  report  for  the  previous 
contract  in  which  this  work  was  initiated. 


LASER- DOPPLER- ANEMOMETER  ( LDA)  MEASUREMENT  AND 
DATA -ACQUISITION  SYSTEM 

LDA  Measurement 

Efforts  were  directed  both  at  installing  the  LDA  apparatus 
in  conjunction  with  the  compressor  test-rig  and  at  developing 
the  LDA  technique  itself. 

LDA  Mounting_on  Compressor_Test  Ri^ 

The  compressor  rig  was  modified  and  added  to  in  order  to 
allow  installation  of  the  LDA  for  velocity  measurements. 

To  allow  optical  access  to  the  flow  passages  of  the 
compressor,  a  2"  x  3  1/2"  window  was  designed  and  installed  on 
the  compressor  rig.  The  final  window  design  involved  l/16"-thick 
Plexiglas  material  cemented  to  a  window  holder.  The  window  holder 
formed  and  held  the  window  material  in  a  cylindrical  shape  matching 
that  of  the  outside  diameter  of  the  compressor-rig's  annular  flow 
passage.  The  window  holder  mates  with  one  of  the  compressor 
rig's  housing  rings.  The  holder  can  be  adjusted  relative  to  the 
housing  ring,  allowing  the  window  to  be  flush  with  the  outside 
diameter  of  the  annular  flow  passage.  The  housing  ring  containing 
the  window  can  be  assembled  into  the  complete  compressor  rig  at 
a  variety  of  axial  stations.  This  allows  LDA  access  not  only  to 
the  compressor  rotor  and  blades  but  also  to  portions  of  the  annular 
flow  passage  at  positions  either  upstream  or  downstream  of  the  rotor 
station.  Optical  calculations  and  tests  with  the  LDA  verified 
that  the  window  would  not  significantly  alter  the  optical  path 


of  the  LDA  beams.  The  possibility  of  LDA  attenuation  due  to 
excessive  reflection  from  the  window  surfaces  was  also  considered. 

It  was  verified  that  commercial  optical  coatings  are  available 
to  minimize  such  reflection,  if  necessary. 

To  provide  support  and  adjustment  of  the  LDA  itself,  a 
table  was  designed  and  constructed.  The  table  locates  the  LDA 
apparatus  approximately  parallel  to  the  compressor  axis, and  an 
adjustable  mirror  provides  redirection  of  the  LDA's  optical  axis 
to  be  transverse  to  the  compressor  axis.  This  arrangement  minimizes 
the  possibility  of  the  LDA  optics  being  struck  in  the  event  of 
a  compressor  blade  failure.  The  design  of  the  LDA  table  provides 
for  three-dimensional  adjustment  of  the  entire  LDA  assembly. 

While  height  adjustment  is  provided,  it  is  not  intended  to  be 
used  as  frequently  as  LDA  position  adjustment  in  a  horizontal 
plane.  Horizontal  adjustment  allows  two  degrees-of- freedom  by 
way  of  an  adapted  milling  table.  The  axes  of  this  table  are 
aligned  to  provide  independent  positioning  of  the  LDA  in  either 
axial  or  transverse  directions  relative  to  the  compressor  rig. 
Sufficient  motion  is  provided  to  allow  use  of  the  entire  window 
in  the  compressor-rig  housing  (approx.  2",  axial)  and  to  locate 
the  LDA  measuring  volume  at  any  radial  position  in  the  annular 
flow  passage  (12"  i.d.,  20"  o.d.).  To  control  vibrations,  the 
LDA  table  was  mounted  on  concrete  pads  keyed  to  the  concrete 
floor  of  the  laboratory.  Initial  qualification  tests  for 
vibration  suggested  that  the  LDA  table  was  not  subject  to  signifi¬ 
cant  vibration  when  the  compressor  rig  was  running. 


Pc.  -.icle_Seeding 


Development  of  the  LDA  technique  primarily  involved  work 
concerned  with  the  particle-seeding  of  the  air  flow  which  is 
required  by  the  technique.  Also  involved  was  development  of 
experimental  methods  and  equipment  to  allow  routine  calibration 
and  performance  checks  of  the  LDA  system. 

The  requisite  characteristics  of  LDA  particles  and  the 
method  of  introducing  them  for  use  was  clarified  by  calculations. 
These  involved: 

(1)  modeling  response  of  the  particle  to  unsteady 
velocities  in  the  surrounding  air, 

(2)  acceleration/deceleration  of  a  particle  injected 
into  a  moving  air  stream, 

(3)  dispersion/dilution  of  particles  during  their 
equilibrium  with  the  surrounding  air  flow,  and 

(4)  concern  for  the  rate  of  data  acquisition  from  LDA 
particles. 


(1)  LDA  Response  to  Unsteady  Air  Velocities 

For  the  purposes  of  estimating  response  to  unsteady  air 
velocities,  the  LDA  particle  was  modeled  as  an  isolated  sphere 
subject  to  Stokes-Law  drag  by  a  surrounding  air  stream.  It  is 
then  possible  to  estimate  the  amplitude  and  phase  error  in 
particle  response.  Under  these  conditions,  the  particle  response 
is  first-order  in  velocity  with  respect  to  time.  Response  can, 
therefore,  be  characterized  in  terms  of  a  time  constant: 


[mV*] 


t  =  1.25  x  10 _5s  =  10  5s 

For  such  a  value  of  the  particle  time  constant,  the  relative 
(velocity)  amplitude  response,  A,  and  relative  phase,  0,  of  the 
particle's  oscillation  can  be  readily  expressed  in  first-order 
form : 

A  =  {1/11  +  ( 2ir f  x )  2  )  }1/2 
0  =  arctan  (2irfT) 

Amplitude  and  phase  characteristics  of  the  baseline  particle  are 
shown  in  Figure  l(for  t  =  10  s)  along  with  additional  curves. 
Response  curves  for  particles  twice  or  half  as  dense  as  water 
(x  =  2  x  10  5s,  t  =  0.5  x  10  5s;  respectively)  and  for 
particles  of  2  urn  and  0.5  ym  diamters  (t  =  4  x  10-5s,  t  =  0.25 
x  10  5s;  respectively)  are  included.  Similar  results,  in 


non-dimensional  form  are  reported  in  Ref.  2  and  are  in  agreement 
with  the  values  of  Figure  1. 


PHASE 


FREQUENCY  (HZ) 


FIGURE  1:  Relative  Amplitude  and  Phase  Lag  (in  degrees) 
Calculated  for  a  Single  Particle  in  Airstream 
of  Sinusoidally  Varying  Velocity 


One  concludes  from  the  values  of  Figure  1  that  liquid 
particle  sizes  near  1  ym  could  be  expected  to  respond  with 
sufficiently  small  error  in  unsteady  velocity  up  to  frequencies 
of  say,  1000  Hz  (A  =  0.997,  0  =  4.3  degrees).  Unsteady  velocity 
components  above  such  a  frequency  are  not  likely  to  be  of  major 
significance  in  the  subject  propagating-stall  situations.  There¬ 
fore,  such  particles  are  suitable,  at  least  on  an  unsteady- 
response  basis,  for  LDA  seeding. 

(2)  Seed-Particle  Equilibration  Following  Injection  Into 
Moving  Air  Stream 

The  need  to  add  or  "seed"  particles  into  the  air  flow  requires 
concern  for  the  speed  with  which  the  seeded  particle  equilibrate 
in  velocity  with  the  surrounding  air  flow.  Assuming  that  Stokes 
Law  of  viscous  drag  applies  (Reynolds  Number  based  on  relative 
velocity  <  0.5),  the  relative  particle  velocity  is  first-order 
in  time  for  a  single,  isolated  particle  suddenly  seeded  (injected) 
into  a  surrounding  air  stream. 

For  the  "baseline"  situation  mentioned  above  (1.  ym  water 
particle  in  STP  air) ,  the  distance  traveled  by  such  a  particle 
in  equilbrating  in  velocity  within  1%  can  be  calculated  (Ref.  2) 
to  be  only  about  3  mm.  On  this  basis,  particle  acceleration/ 
deceleration  allows  that  particle  injection  can  be  at  virtually 
any  practical  distance  upstream  of  the  LDA  measuring  volume.  That 
is,  factors  other  than  simple  particle  acceleration/deceleration 
will  dominate  choice  of  an  appropriate  upstream  position  for 
seeding. 


(3)  Carrier-Gas  Equilibrium  Following  Injection 

LDA-particle  seeding  involves  injection  of  particles  and, 
with  them,  a  surrounding  carrier-gas  (air)  flow.  For  this  reason, 
an  additional  aspect  of  particle  response  involves  equilibration 
of  the  carrier-gas  flow  with  the  main-stream  air  flow.  This 
situation  was  modelled  as  a  turbulent- jet-mixing  problem  in  order 
to  estimate  the  distance  of  particle  (jet)  travel  required  between 
the  injection  point  and  measuring  point  to  insure  near-equilibration 
of  velocities. 

Jet-mixing  theory  by  Abramovich  (Ref.  3)  was  applied  to 
the  problem  of  calculating  the  distance  downstream  of  carrier-gas 
injection  at  which  the  carrier-gas  (and  particle)  velocity  and 
the  free-stream  velocity  would  be  equilibrated  within  a  specified 
difference.  While  the  unsteady  velocities  expected  in  turbomachine 
tests  are  probably  not  well  approximated  by  such  a  jet-mixing 
model,  the  intent  of  the  modelling  was  to  establish  whether 
carrier-gas  equilibration  was  likely  to  represent  a  practical 
problem  or  not;  the  jet-mixing  approach  was  thought  to  provide 
information  as  to  upper  limits  on  the  distances  downstream 
required  for  equilibration  and  a  basis  for  judging  whether  this 
problem  might  be  practically  significant  or  not. 

The  results  of  calculations  showed  that  unless  the  carrier- 
gas  and  (mean) free-stream  velocities  are  well  matched  at  the  point 
of  carrier-gas  (and  particle)  injection,  large  distances  downstream 
are  required  (e..g,  1000  nozzle  diameters)  and  considerable  dilution 
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and  jet  spreading  (e.g.,  to  tens  of  nozzle  diameters)  of  the 
carrier-gas  will  have  occurred  by  the  point  where  the  carrier-gas 
and  free-stream  have  nearly  equilibrated. 

It  was  concluded  that  particular  attention  would  have  to 
be  paid  to  this  factor  during  experimental  development  of  the  seed¬ 
ing  technique. 

(4)  LDA  Data  Rate 

Particularly  in  light  of  a  possible  data-reduction  scheme 
requiring  high  sampling  rates  (Ref.  1) ,  possible  LDA  data  rates 
were  estimated  in  order  to  check  feasibility. 

Theoretically,  data  rate  is  limited  by  the  fringe  spacing 
in  the  LDA  measuring  volume,  the  counting  algorithm  employed  by 
the  LDA  processor,  the  seed-particle  speed  (air  speed) ,  and  the 
seed-particle  concentration.  Practically,  the  nature  of  the 
particle  (size,  shape,  material)  and  the  laser-beam  intensity 
also  influence  actual  data  rates. 

Theoretical  limits  on  data  rates  are  shown  in  Figure  2 
where  fringe  spacing  and  particle  velocity  are  shown  as  determining 
maximum  data  rate.  The  data  presented  are  based  on  an  assumed 
eight  fringe  crossings  being  required  for  a  single  measurement. 

It  is  clear  from  the  results  that  a  reasonable  velocity  range 
may  be  covered  with  data  rates  high  enough  to  allow  gated  LDA 
measurements  of  the  sort  under  consideration  for  the  present  work. 


2eXeA°Em®nt.s_to  Routinize  LDA  Use 


Many  bench-top  experiments  were  carried  out  in  the  course 
of  developing  the  LDA  methods.  A  major  portion  of  these  were 
aimed  at  providing  effective  seeding,  flow  situations,  etc.  for 
convenient,  effective  checks  on  operation,  calibration,  etc.  of 
the  LDA.  The  earliest  of  these  experiments  typically  involved 
air  flows  and  smoke  particles,  in  anticipation  of  the  need  for  a 
somewhat  similar  method  for  use  with  the  compressor  rig.  These 
were,  however,  abandoned  in  favor  of  latex  spheres  in  water  and 
a  water  "tunnel"  to  provide  more  convenient  reproducibility. 

LDA  Data  Acquisition  &  Storage 

Both  hardware  development  and  software  programming  were 
undertaken  during  the  subject  contract  period  with  the  aim  of 
providing  for  gated  signals  from  the  LDA.  Figure  3  shows 
schematically  the  anticipated  "sparse-burst"  nature  of  the  LDA 
data  to  be  acquired  and  stored.  The  data  bursts  shown  occur 
during  the  short  "gate-open"  periods  for  which  a  selected  portion 
of  a  blade  passage  is  within  the  measuring  volume  of  the  LDA. 

LDA  velocity  measurements  obtained  by  the  LDA  at  other  times  are 
not  of  interest  and  therefore  are  blanked  out  during  the  relatively 
long  gate-closed  periods. 

To  provide  for  the  expected  data,  the  data-acquisition  and 
storage  system  shown  schematically  in  Figure  4  was  designed  and 
construction  initiated. 
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The  backplane  wiring  of  a  card  rack  for  the  single-board 
computer  and  two  additional  boards  was  completed  (with  checking) . 
One  of  the  additional  boards,  containing  clock  and  DMA  controller 
was  constructed.  However,  at  completion  of  the  contract,  the 
second  board  containing  the  serial  interface  had  not  yet  been 
contructed. 

Programming  of  the  microcomputer  to  provide  for  gating  of 
the  LDA  signal,  etc.  was  started  using  a  microprocessor  development 
system,  but  this  work  was  not  completed. 

Figure  5  shows  in  more  detail  the  functional  blocks  which 
comprise  the  system.  To  accommodate  the  high  data  rates  desired, 
direct  memory  access  (DMA)  was  provided  under  the  control  of  a 
single-board  microcomputer  based  on  an  Intel  8085A  microprocessor. 
The  microprocessor  receives  once  per  revolution  pulses  from  a 
sensor/signal  conditioner  installed  on  the  compressor  test  rig. 
These  pulses  serve  as  the  basis  for  the  microcomputer  to  determine 
at  what  time  the  gate  should  be  opened,  allowing  any  data  acquired 
from  the  LDA  to  be  transmitted  and  stored  in  memory.  The  micro¬ 
processor,  in  turn,  signals  the  DMA  controller  at  the  appropriate 
time,  and  data  from  the  LDA  Counter  Processor  (fD,  the  Doppler 
frequency,  corresponding  with  measured  velocity)  begins  to 
flow  into  the  memory  of  the  system.  For  each  12-bit  item  of 
velocity  data  from  the  LDA  that  is  stored,  a  corresponding  12-bit, 
100-KHz  clock  reading  is  stored  to  record  the  time  (t)  at  which 
the  velocity  was  measured.  The  clock  is  effectively  read  with 
24-bit  resolution,  however,  by  storing  either  most-  or  least- 
significant  sets  of  12  bits.  Most-significant-bit  (MSB)  data  are 
stored  only  when  the  MSB  data  are  incremented.  This  is  signaled 


fD  •  DOPPLER  FREQUENCY 


t  -  TIME 


FIGURE  5:  Functional  Block  Diagram  of  Data  Acquisition 
and  Storage  System 
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by  a  code  value  appearing  in  the  corresponding  Doppler  frequency 


register  rather  than  an  actual  vlaue  of  fp.  The  12-bit  data  are 
stored  in  banks  of  RAM  for  later  recall  and  processing.  Since 
both  LDA  output  and  timing  clock  data  are  in  coded  form,  some 
processing  is  required  merely  to  retrieve  the  physically  meaningful 
values  from  the  stored  data. 

The  microcomputer  serves  to  withdraw  data  from  memory  via 
the  DMA  controller  and  transmit  the  data  to  another  data-processing 
computer.  This  transmission  is  designed  to  be  via  an  RS-232 
serial  port. 


LDA  DATA  PROCESSING 

Numerical  simulations  were  continued  with  the  aim  of 
evaluating  alternative  schemes  for  deducing  unsteady  velocities, 
power  spectra,  etc.  from  the  previously  described  LDA  data. 

The  work  previously  reported  (Ref.  1)  was  extended  to  determine 
the  feasibility  of  signal  reconstruction  for  the  present  application 
according  to  the  method  of  Yen  (Ref.  4).  Interest  in  this  method 
arose  from  earlier  simulations  by  other  methods  (Ref.  1)  which 
indicated  numerical  difficulties  with  reconstructing  the  original 
signal  from  the  sparse-burst  samples  of  the  subject  experiment. 

Simulations  using  Yen's  Minimum  Energy  Method 

The  "minimum  energy"  method  of  Yen  (Ref.  4)  was  tested 
with  the  hope  of  finding  that  it  would  minimize  the  tendency 
observed  with  other  methods  to  generate  large,  false  peaks  or  "blips 
in  the  reconstructed  signal.  These  spikes  presumably  derive  from 


numerical  inaccuracies. 


Yen's  method  involves  numerically  determining  a  fictitious 
set  of  sample  values,  equally  spaced  in  time  that  would,  upon  signal 
reconstruction,  match  the  actual  measured  signal  values  at  those 
arbitrary  times  at  which  the  actual  samples  were  taken.  The 
matching  is  subjected  to  the  constraint  that  the  sum  of  the 
squares  of  the  sample  values  of  the  fictitious  set  should  be  a 
minimum. 

Figure  6  shows  an  example  of  a  mixed  frequency  signal  from 
which  samples  were  drawn  in  bursts  (see  Fig.  3) ,  with  two  samples 
randomly  timed  within  each  burst  or  gate-open  period.  The  points 
shown  represent  the  pairs  of  samples.  Yen's  algorithm  was  then 
used  to  reconstruct  the  signal  from  these  pairs  of  samples  by 
using  the  intermediary  of  the  fictitious,  uniformly  spaced  series 
of  samples  referred  to  above. 

Figure  7  shows  the  signal  reconstructed  from  the  pairs  of 
randomly  spaced  samples  for  an  example  case  of  80  samples  (40 
pairs) .  At  two  times  (approx.  29  and  39  time  units)  large  spikes 
can  be  seen,  and  in  fact,  smaller  spikes  are  apparent  in  the 
neighborhood  of  every  pair  of  samples.  All  such  spikes  are 
apparently  the  result  of  numerical  errors,  probably  deriving  from 
inaccuracies  in  the  matrix  inversion  algorithm  used  to  solve  the 
reconstruction  equations.  This  is  evidenced  by  the  fact  that  the 
occurrence  of  these  spikes  and  their  magnitudes  differed  considerably 
from  one  matrix- inversion  algorithm  to  another.  For  example, 
an  early  Scientific  Subroutine  Package  (SSP)  program  was  found 
to  give  somewhat  less  trouble  than  several  more-modern  International 
Mathematics  and  Statistics  Library  (IMSL)  routines  for  matrix  inversion 


r lgina 


It  appears  that  low-pass  filtering  of  the  reconstructed 
signal  could  readily  eliminate  these  spikes  without  interfering 
appreciably  with  the  much- lower- frequency  components  which 
represent  the  original  signal.  Figures  8a  and  8b  show  power 
spectra  for  the  original  and  the  reconstructed  signals  of  Figures 
6  and  7.  It  is  clear  that  the  obvious  high-frequency  spikes  of 
the  reconstructed  signal  do  not  seriously  impede  meaningful 
spectral  analysis  at  the  lower  frequencies  of  interest.  However, 
it  is  notable  that  in  all  cases  tested,  increasing  the  number  of 
samples  always  led  to  increasingly  large  spikes  occuring  at  larger 
sample  sizes.  Thus,  there  is  evidence  of  an  inherent  mismatch 
between  the  reconstruction  method  and  the  drive  for  increased 
resolution  and  signal-to-noise  ration  via  increased  sample  size. 

In  cases  involving  a  relatively  small  number  o!;  samples 
(e.g.,  10),  spikes  were  not  observable  in  reconstructed  signals 
and  certainly  were  not  troublesome.  However,  either  increasing 
the  number  of  samples  in  each  burst  (during  each  gate-open  period) 
or  increasing  the  number  of  gate-open  periods  led  to  increasingly 
large  spikes,  finally  to  the  point  of  seriously  disrupting  the 
reconstructed  signal. 

The  spikes  observed  in  the  reconstructed  singals  were 
studied  as  to  their  sensitivity  to  the  parameters  defining  the 
sampling  situation,  especially  the  time  duration  of  the  "gate-open" 
period  relative  to  the  duration  betwen  gate-open  periods.  At 
gate-open  periods  sufficiently  large  relative  to  the  gating  period 
(e.g.  30%),  it  was  found  that  the  spikes  virtually  disappeared. 


In  light  of  the  numerical  difficulties  encountered  with 
implementing  Yen's  method  of  minimum  energy,  the  typical  form 
of  the  matrices  involved  was  studied.  The  algorithm  for  re¬ 
construction  involves  inversion  of  a  matrix  with  elements  near 
unity  in  value  on  the  diagonal  and  near  the  diagonal.  Elements 
of  the  matrix  are  much  smaller  elsewhere  (far  off-diagonal) . 

This  approximately  "block  diagonal"  form  with  near-zero  elements 
off-diagonal  would  tend  toward  a  singular  matrix  of  zero  determinant 
in  the  limit  of  very  narrow  gate-open  periods  with  much  longer 
intervening  gate-closed  periods.  It  is  possible  that  specialized 
inversion  methods  for  such  block  matrices  are  available.  However,  it  is 
clear  that  the  near-singularity  of  these  matrices,  worsening  as 
gate-open  times  decrease,  represents  a  major  practical  difficulty 
in  implementing  Yen's  method. 


Alternative  Algorithms 

In  light  of  the  difficulty  in  matrix  inversion  which 
seems  inherent  in  Yen's  method,  it  is  worthwhile  to  consider 
algorithms  which  are  better  adapted  to  the  sampling  required  by 
the  present  work,  i.e.,  very  short  gate-open  periods.  Yen's 
method  apparently  works  best  in  the  limit  of  long  gate-open 
periods  (approaching  100%) ,  and  this  limit  is  the  opposite  of 
that  toward  which  it  is  desirable  to  strive  in  the  present  case. 

It  is  potentially  advantageous  to  adopt  an  algorithm  which  shows 
promise  of  applicability  in  the  opposite  limit,  that  of  very  short 
gate-open  periods. 


An  example  of  a  procedure  which  is  somewhat  better  adapted 
than  Yen's  to  the  present  situation  is  that  of  using  pairs  of 
sample  data  within  gate-open  windows.  From  these,  mean  and  first- 
derivative  values  may  be  calculated.  An  existing  explicit  algorithm 
(Ref.  5)  is  then  feasible  for  reconstructing  signals  from  such 
signal  values  and  first  derivatives  which  are  uniformly  spaced 
in  time. 

Alternatively,  pairs  of  sample  data  might  be  used  to 
calculate  means  and  first  derivatives  to  which  fictitious 
uniformly  spaced  "samples"  could  be  fitted  and  processed. 

Again,  for  these  pairs,  an  explicit  algorithm  which  requires 
uniformly  spaced  rather  than  randomly  spaced  samples  within 
each  gate-open  period  is  available  (Ref.  5) . 

Unfortunately,  the  only  evident  alternatives  are  not 
likely  to  be  useful  if  desirably  high  sampling  rates  are  achieved. 

If  more  than  a  pair  or  perhaps  three  samples  are  acquired  during 
single  gate-open  period,  these  explicit  methods  are  apparently 
not  applicable  and  no  others  are  known.  Still,  it  is  just  such 
high  rates  that  would  provide  the  greatest  potential  for  determining 
unsteady  velocity  components  at  substantial  multiples  of  the 
rotational  frequency  of  the  rotor,  i.e.,  for  high  engine  orders. 
Multiple  samples  per  gate-open  period  increase  the  Nyquist 
frequency  for  the  sampling  to  an  equal  multiple  of  one-half  the 
rotational  frequency.  More  than  two  or  three  samples  per  gate-open 
period  are  necessary  in  order  to  deal  with  unsteady  components 


much  above  the  rotational  frequency.  Extension  of  the  gated-LDA 
approach  to  measuring  unsteady  frequency  components  at  the  relative 
ly  high  frequencies  involved  in  multi-patch  stall  phenomena  seems, 
therefore,  unlikely. 


ANALYTICAL/ CFD  MODEL  OF  CASCADES 


Coupled  with  the  experimental  portion  of  the  research  there 


has  been  an  effort  to  construct  a  tentative  analytical  model  of 


the  time-varying  behavior  of  a  cascade  of  airfoils.  The  approach 


has  been  to  treat  a  cascade  as  a  system  of  N  interfering  airfoils 


without  any  a  priori  assumptions  concerning  uniformity  of  loading 


from  airfoil  to  airfoil.  A  search  for  other-than-unif orm-f low 


solutions  to  the  model  equations  should  yield  some  insight  into 


flow  about  highly  loaded  cascades. 


Arbitrarv  Loading  Cascade  Model 


The  infinite  two-dimensional  cascade,  obtained  by  unwrapping 


an  annular  cascade  of  differential  radial  extent,  will  have  periodi¬ 


city  imposed  over  every  N-bladed  group  of  airfoils  (Fig.  9  ) .  Here 


N  represents  the  number  of  blades  in  the  annular  cascade  from 


which  the  two-dimensional  cascade  was  developed.  Typically  10  <  N  <  100 


and  is  often  chosen  to  be  prime.  For  simplicity  the  airfoils  are 


uncambered*  and  thin  so  that  only  the  velocities  normal  to  the 


chords  need  be  considered  in  the  boundary  condition. 


Under  these  proscriptions,  including  incompressibility  of  the 


fluid,  the  statement  of  the  boundary  conditions  leads  to  a  matrix 


expression  of  the  following  form 
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ightly  loaded,  slightly  cambered  airfoils  the  idemvector  on  the 
RHS  of  Eq.  (1)  would  be  augmented  by  a  function  of  x. 
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FIGURE  9:  Geometry  of  infinite  cascade  derived  from  N-bladed  annular 
cascade. 


where  Vmn  are  integral  operators,  Vnn  are  singular  integral  operators 

(i.e.  on  the  main  diagonal)  and  Yn(S)  in  the  unknown  vorticity,  or 

distributed  circulation,  on  the  nth  airfoil,  1  <  n  <  N.  V  is  the 

—  —  m 

component  of  the  mean  velocity  normal  to  chordlines  and  is  determined 
by  specifying  the  (constant)  vector  velocity  at  upstream  infinity 
and  subtracting  from  it  the  velocity  induced  by  all  the  vortices 
(both  free  and  bound)  in  the  field.  In  addition  there  are  free 
vortices  in  the  field.  The  ith  free  vortex  of  strength  r.  and 
located  at  and  its  N-spaced  images  induce  the  velocity,  normal 

to  the  chord  of  the  nth  airfoil  given  by  r\(£^,  n^)vn^(x). 

The  integral  operators  ,  accounting  as  they  do  for  induction 
effects  due  to  every  Nth  airfoil  in  an  infinite  array,  are  given 
by  (see  Fig.  for  the  geometry  of  the  cascade  expressed  in  non- 


orthogonal  coordinates) . 


-)  Kmn(x  -  £)  d£ 


where 


K  =  ~  i  _ £-x  +  (n-rrH-qN)s' _ 

mn  2V  q=-°°  [E-x  +  (n-m+qN)s']2  +  [n-m+qN)s"]2 

s'  =  s  sin  B,  s"  =  s  cos  S 

Thus,  V  „  operates  on  y  to  produce  a  contribution  to  the 
mn  n 

velocity  at  the  mth  blade. 

Similarly,  the  velocities  induced  by  the  discrete  vortices 


are  given  by 


v  .  =  —  I 
mi  2  tt 


-  x  +  sinB  -  ms'  +  qNs' 


2 

[£^-x  +  ni  sinB  -  (m-qN)s')  +  [n^  cos  8  -  ms"  +  qNs 


Since  there  seems  to  be  no  direct  method  of  manipulating  the 


matrix  of  operators,  the  equations  are  solved  by  assigning  an 
expansion  of  aerodynamic  modes  to  each  vorticity  distribution,  say 

-  Vir!^  +  +  Vi-t/b)1*  a+c/b)*5  < 

+  E  S/bd-S/b)15  [l+Z/b)h,  -b  <£  <  b 
n  —  — 

then  each  Nth  airfoil  will  have  the  same  distribution,  i.e. 

Yn  =  Yn+N  =  Yn-N  =  Yn+2N  etc‘ 

The  coefficients  An  and  Bn  allow  for  singularities  at  the  leading 
and  trailing  edges,  respectively  whilst  Dn  allows  for  an  adjustment 
of  the  circulation.  En  provides  for  biasing  the  distribution  more 
to  the  leading  or  trailing  edges,  depending  on  sign. 

Then  each  integral  is  evaluated  at  a  number  of  points,  x,  along 
the  chord  of  each  airfoil:  four  points  in  the  case  of  the  example 
expansion  in  Eq.  (6) .  This  is  tantamount  to  satisfying  the  boundary 
conditions  by  collocation,  so  that  with  the  current  choice  of  4 
terms  in  Eq.  (6),  values  of  x  =  -b,  -b/2,  b/2,  b  would  be  appropriate 

With  these  definitions  the  equations  are  reduced  to  true 
matrix  equations  of  order  4Nx4N  (in  this  particular  example). 

[a]  {A}  =  -V  {1}  -IT.  {v.}  (7 

m  ,ii 

l 

Here  {A}  is  the  column  vector  of  unknown  vorticity  coefficients 

{A}  =  [A^,  B^,  D^,  E^,  A2,  B2,  D2,  A3 ,  .  .  .  ,  D^]  (8 

and  Vm  is  a  constant  whose  physical  meaning  has  been  described. 

In  addition  there  may  be  a  number  of  free  vortices,  in  the  field, 
and  their  N-spaced  images  (see  Fig.  1)  which  also  induce  the 
velocities  expressed  by 


Fi  ^i' ni  [ vi i '  vli'1)'  vli(2)'  vii'3)'  v2i(0),  v2i(1)' 


.vNi(3)] 


Here  are  the  coordinates  of  the  vortex  of  strength  F^. 

The  coefficient  matrix  [a]  is  of  order  4NX4N  and  has  the 
following  general  form 
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The  superscripts  0,1,2  and  3  refer  to  the  collocation  points  chosen 
for  satisfying  the  boundary  conditions. 

The  coefficient  matrix  is  the  result  of  inserting  Eqs.  (6)  and 
(3)  into  Eq.  (2)  and  performing  the  indicated  operation,  e  .g  . 

«o  .  ,b  .bis.i,  J_  •  - i-  <2-uWs‘ - 

12  _h‘  'b-v  2ri  q,.„[£-(-b/2)  +  <2-l+qN)s'l  +  [ (2-1+qN) s" ] * 


Referring  to  Eq.  (5),  the  velocity  induced  at  x  =  -b/2  on  airfoil 
n  -  1  by  the  ith  free  vortex  is  given  by 

n  1  00  -(-b/2)  +  nA  sine  -  l.s'  +  qNs' 


(-b/2)  +  n. sing  -  (m-qN)s']  +  [n.cosg  -(m-qN)s"] 


Nonuni  form  distribution  which  result  from  solving  Eq.  (7)  for 

the  vector  {A}  will  depend  critically  on  the  field  of  free 

vortices  I { r . } . 
i 


Uniform  Cascades  as  a  Special  Case 

In  the  absence  of  free  vortices,  {I\}  =  {0},  the  solution 
of  Eq.  (7)  will  yield  the  typical  result  of  uniform  loading  if 
there  is  imposed  the  same  Kutta  condition  of  smooth  flow  off  the 
trailing  edge  of  each  blade.  Setting  all  Bn  =  0,  so  that  the 
distributed  circulation  drops  to  zero  at  the  trailing  edge  of 
every  blade,  results  in  more  simultaneous  equation  than  unknowns. 

Every  4th  column,  starting  with  the  2nd,  (i.e.,  columns  2,6,10 
etc.)  may  be  eliminated  since  these  terms  always  multiply  a 
Bn  component.  Furthermore,  the  number  of  rows  may  also  be  reduced 
by  one  quarter  by  eliminating  one  of  the  collocation  points,  say 
x  =  -b,  the  leading  edges  of  the  foils. 

It  should  be  emphasized  that  these  consequences  of  applying 
the  Kutta  condition  are  without  mathematical  foundation.  This 
"condition"  is  an  observable  physical  phenomenon  which  can  be 
explained  in  the  light  of  the  slight,  but  nonvanishing,  viscosity 
of  the  fluid.  Thus,  in  any  inviscid  approximation  to  the  flow  such 
as  the  present  one  the  Kutta  condition  becomes  an  arbitrary  mathematical 
statement.  The  lesson  to  be  learned  from  the  consideration  is 
that  dynamic  stalling  of  the  cascaded  airfoils  cannot  be  determined 
entirely  without  the  specification  of  analogous  "Kutta  conditions". 

When  the  passage  flow  reverses,  smooth  flow  off  the  leading  edges 
may  be  expected  at  certain  times.  If  temporally  separated  flow 


appears  with  separation  streamline,  emanating  from  leading  and/or  trail 
ing  edges,  additional  "separation  conditions"  may  have  to  be  imposed. 

The  partitioning  of  [a]  from  Eq.  (9)  into  submatrices, 

of  order  3x3,  reveals  certain  symmetries  when  each  submatrix 

is  treated  as  an  element,  a 

'  mn 

i)  All  elements  on  the  main  diagonal  are  equal,  i.e. 
all  ann  are  equal  irrespective  of  n. 

ii)  All  elements  in  a  line  parallel  to  the  main  diagonal 
are  equal  to  each  other,  i.e.  depends  only  on  m-n 

and  not  on  m  and  n  individually.  (Statement  ii  includes  i) . 

iii)  If  m  and  n  for  one  element  are  such  that  their  difference 
is  N  greater  or  less  than  m-n  for  another  element,  the 
two  terms  are  equal. 

All  these  properties  taken  together  are  such  that  the  matrix  of 
subelements  is  a  circulant  matrix. 

Provided  that  [a]  is  non-singular,  this  results  in  a  solution 
to  Eq.  (1)  in  which  all  the  An,  Dn ,  En  triplets  are  equal.  The 
loading  is  the  same  on  each  airfoil.  Since  the  loading  is  found 
to  be  uniform  in  this  case,  the  a  priori  assumption  that  all  the 
free  vortices,  I\,  are  absent  has  been  justified. 

Nonuniform  Solutions 

Clearly,  a  preliminary  estimate  of  nonuniform  loading  flow 

solutions  may  be  synthesized  by  setting  An  =  0  for  certain  values 

of  n  and  B  =  0  for  other  values  of  n.  Simultaneously  it  is 
n 

necessary  to  eliminate  the  column  in  [a]  associated  with  each 


null  coefficient.  For  every  column  eliminated  it  is  also  necessary 
to  eliminate  a  row  in  order  to  maintain  [a]  as  a  square  matrix. 

The  resulting  stall  "patches"  might  then  be  expected  to  appear 
in  the  vicinity  of  those  blades  having  An  =  0. 

This  solution  is  somewhat  unrealistic  in  the  sense  that 
there  are  no  free  vortices  downstream  of  the  cascade,  nor  do  the 
stall  patches  propagate.  These  observations  are  tantamount  to 
stating  that  the  flow  is  steady  in  contradiction  to  what  should 
be  the  physical  reality. 

More  realistic  solutions  of  the  equations  will  depend  upon 
the  ability  to  link  the  temporal  fluctuation  in  the  An(t) ,  Bn(t) 
and  Dn(t)  to  the  strengths  of  the  free  vortices  I\  and  their 
trajectories.  This  is  evident  from  the  fact  that  the  circulation 
around  the  nth  airfoil  is  given  by 

b 

Circn  =  /  y(0  d£  =  ub(An  +  Bn  -  1/2  Dn)  (13 

As  the  coefficients  A  ,  B  and  D  change  with  time,  the  circulation 
changes  and  the  shed  circulation  appears  in  the  form  of  the  free 
vortices  in  order  to  conserve  the  total  circulation.  Since  the 
trajectory  of  each  free  vortex  depends  on  time  and  on  the  flow 
velocity  at  the  location  of  the  vortex,  it  is  clear  that  one  method 
of  obtaining  the  instantaneous  field  of  velocities  (and  locus  of 
vortices)  is  by  a  time-marching  algorithm  starting  with  some  assumed 
initial  configuration. 

Although  there  are  some  avenues  to  be  investigated  in  attempting 
to  obtain  the  more  realistic  solutions,  these  have  not  been  pursued 
up  to  the  time  of  concluding  this  contract. 


RESEARCH  COMPRESSOR 


This  apparatus,  constructed  with  ONR/NavAir  funding, 
is  a  single  stage  axial  flow  compressor  embodying  a  number 
of  distinctive  features.  Some  of  these  features  are  unique: 

.  Constant  flow  annulus  has  a  20  inch  (0.508m)  tip  diameter  and 
12  inch  (0.305m)  hub  diameter. 

.  Rotor  drum  is  fitted  with  36  uniform  (chord  5.08cm  and  aspect 
ratio  of  3)  untwisted  blades.  Blades  are  shaft-mounted  so 
that  stagger  may  be  changed  by  resetting  a  unison  ring  within 
the  drum. 

.  Key  housing  sections  are  fabricated  from  polymethyl  methacrylate 
(PMM)  for  ease  of  flow  visualization.  Housing  sections  de¬ 
signed  for  accepting  LDA  laser  beams  are  of  cast  aluminimum 
with  small  PMM  windows  set  in  flush  with  the  housing  I.D. 

.  Inlet  guide  vanes  and  exit  guide  vanes  (stator  blades) are 

variable  by  means  of  unison  rings  and  at  design  theoretically 
provide  constant  incidence  on  the  rotor  blades. 

.  Discharge  throttle  is  followed  by  a  large  shipboard  ventilation 
blower  to  allow  adjustment  of  throughflow  independently  of  com¬ 
pressor  speed. 

.  Variable  speed  drive  is  provided  by  a  20  h.p.  (14.9  kW)  D.C. 
motor  with  solid  state  controller  and  feedback  tachometer 
driving  through  a  belted  pulley  system. 

.  Intake  is  conventional  bellmouth  and  nose  bullet  with  four 
front  main  bearing  support  struts. 

The  compressor  has  been  entirely  assembled  and  run  at  a 
number  of  speeds  to  calibrate  the  controller/tachometer.  Rotor 
blade  setting  was  at  45  degrees,  exit  guide  waves  were  at  design. 


Since  delays  with  the  laboratory  program  had  prevented  the  pro¬ 
curement  of  any  LDA  data  on  the  compressor  directly,  some  rudi¬ 
mentary  measurements  were  made  using  alternative  instrumentation 
in  the  last  few  weeks  of  the  contractual  period. 

At  the  nominal  setting  of  the  inlet  guide  vanes,  the 
midradius  blade  angle  is  zero.  Various  rotor  speeds  and  throttle 
settings  were  used  in  the  first  tests  using  a  wedge  probe  to 
determine  the  velocity  triangles  entering  the  rotor  at  this 
reference  radius.  At  combinations  of  speed  and  throttle  where 
the  rotor  incidence  was  in  excess  of  about  11  deg.  it  was 
found  that  the  flow  was  increasingly  unsteady  and  accompanied  by 
an  audible  change  in  the  acoustic  "signature".  The  presumption 
that  the  propagating  stall  was  present  was  provisionally  confirmed 
with  a  small  linen  thread  tuft  in  the  flow  which  became  increasing 
agitated  as  the  incidence  increased  above  about  11  deg.  (It  was 
the  objective  of  the  entire  research  program  to  replace  such 
provisional  and  qualitative  observations  by  quantitative  in-rotor 
determinations  of  the  unsteady  velocity  field) . 

Additional  tests  with  the  fan  on  or  off,  and  with  different 
settings  of  the  inlet  guide  vanes  (from  -9°  to  15°  from  nominal) 
confirmed  that  rotating  stall  was  present  when  the  incidence,  or 
rotor  blade  loading,  was  increased  above  some  critical  value. 

Two  additional  observations  were  also  possible  in  the  limited  final 
testing  under  the  contract:  (1)  the  region  of  large  incidence  and 
large  flow  fluctuations  (the  stall  patches)  appeared  first  and  were 
most  intense  at  the  rotor  blade  tips,  and  (2)  negative  axial  flow 


CONCLUSIONS 


The  work  leading  up  to  this  final  report  has  been  focussed 
primarily  on  the  empirical  mapping  of  the  unsteady  flow  field  in 
the  rotor  passages  of  an  axial  flow  compressor  (single  stage) . 
Although  this  experimental  result  was  not  finally  achieved  within 
the  contract  time  limit,  several  important  conclusions  were 
reached  in  relation  to  the  data  acquisition  system,  the  data 
processing  algorithms,  the  analytical/CFD  modeling  of  compressor 
stall  and  the  research  compressor  itself. 

Design  of  the  data  acquisition  system  was  governed  by 
the  presumed  necessity  of  following  the  evolution  of  the  stall 
process  within  a  single  (the  same)  blade  passage.  With  the  fixed 
LDV  this  implied  that  velocity  realizations  could  be  obtained  only 
during  a  small  fraction  of  each  rotor  revolution.  Hence  the  data 
sample  on  which  a  spectral  analysis  could  be  performed  was  extremely 
sparse,  the  more  so  since  each  passage  width  had  to  be  subdivided 
into  a  number  of  still  smaller  intervals  (windows) .  In  order  to 
characterize  the  unsteady  velocity  distribution  at  least  ten  and 
probably  20  or  30  subdivisions  will  be  required.  Microprocessor- 
based  hardware  and  software  have  been  developed  that  will  theoreti¬ 
cally  accept  this  quality  of  data  and  output  the  desired  unsteady 
flow  field  information.  Simulations  have  demonstrated  this  con¬ 
clusion. 

Seeding  of  the  compressor  flow  to  increase  the  data  rate 
will  be  required.  Some  studies  reported  herein  have  indicated, 
provided  the  carrier  gas  and  free  stream  velocities  are  fairly 
well-matched,  that  the  seeding  can  be  accomplished  and  will  produce 


the  desired  result.  A  residual  concern  is  the  health  hazard  or 
the  particular  seeding  material  that  is  to  be  used.  A  tentative 
candidate  material  is  glycerine. 

These  preliminary  results  having  to  do  with  the  LDA 
measurements  have  led  to  the  major  conclusion  that  the  experimental 
data,  and  its  spectral  analysis,  should  henceforth  be  based  on 
more  than  one  (preferably  all)  of  the  flow  passages.  Hence  the 
emphasis  should  shift  from  gating,  per  se ,  to  an  accurate  time/space 
location  of  any  and  all  velocity  realizations  that  are  obtained. 

The  microprocessor  data  reduction  algorithm  will  then  be  based  on 
the  presumed  identity, albeit  shifted  in  time,  of  the  processes  in 
all  the  blade  passages. 

The  flow  modelling  initiated  toward  the  end  of  the  contract 
period  was  sufficient  to  demonstrate  that  a  nonuniform  distribution 
of  cascade  loading,  including  local  regions  of  reverse  flow,  is 
possible  and  satisfies  all  the  governing  boundary  conditions. 
Questions  of  flow  stability  and  rate  of  propagation  of  the  non¬ 
uniformity  (stall  patches)  cannot  be  settled  without  augmenting 
the  model.  In  the  future,  it  will  be  necessary  to  add  a  field  of 
free  vorticity,  in  addition  to  the  bound  vorticity  studied  here, 
in  order  to  predict  stability  and  compute  propagation  rates. 

The  method  for  doing  this  will  involve  the  application  of  the 
Random  Vortex  Method,  or  a  similarly-based  algorithm,  to  the  flow 
near  and  downstream  of  the  cascade.  It  should  also  be  possible  to 
qualify  or  validate  the  model  on  a  linear  cascade  of  limited  extent, 
or  possibly  a  stationary  annular  cascade.  However,  the  single-stage 


axial-flow  research  compressor  and  LDA  development  should  continue 
for  ultimate  corroboration  of  the  results  in  a  more  realistic, 
engine-like  geometry. 

The  eventual,  long-term  objective  is  to  be  able  to  predict 
the  nunber  of  stall  zones,  and  their  tangential  rate  of  propagation 
for  application  to  axial  compressors  and  fans.  This  information  is 
critical  for  the  prediction  of  blade  resonance  in  rotating  stall, 
for  the  prediction  of  the  inception  of  stall  flutter  and  possibly 
also  for  the  description  of  compressor  surge.  If  these  goals  can 
be  attained,  they  will  represent  the  first  instances  in  which  the 
phenomena  were  accurately  described  on  the  basis  of  unsteady 
cascade  aerodynamics  rather  than  stage  characteristics. 


REFERENCES 


Sisto,  F. ,  "Research  on  the  Behavior  of  Cascaded  Airfoils 
Under  Conditions  of  High  Mean  Loading  and  Flow  Unsteadiness," 

Final  Technical  Report  on  contract  Noool4-79-C-0765  NR-094-391, 
Stevens  Institute  of  Technology,  30  June  1982. 

Dring,  R.P.,  "Sizing  Criteria  for  Laser  Anemometry  Particles", 

J.  of  Fluids  Engineering  (ASME) ,  Vol.  104,  March  1982,  pp.  15-17. 

Abramovich,  G.N.,  The  Theory  of  Turbulent  Jets,  M.I.T.  Press, 
Cambridge  (1963),  Chapt.  5. 

Yen,  T.L.,  "On  Non-Uniform  Sampling  of  Bandwidth-Limited  Signals", 
IRE  Trans.  Grant  Dreory,  Vol.  CT-3,  Dec.  1956,  pp.  251-257. 

Papoulis,  A.,  Signal  Analysis,  McGraw-Hill,  N.Y.  (1977),  pp.  191-196 


REPORT  DOCUMENTATION  PAGE 


1  REPORT  NUMBER 

ME- RT- 8 3 007 


4  Title  fnd  Submit)  5.  TYPE  OF  REPORT  &  PERlOC  covered 

Behavior  of  Cascaded  Rotor  Airfoils  Under  Final  Technical  Report 
Conditions  of  High  Mean  Loading  and  Flow  for  period  6/1/83-5/31/ 
Unsteadiness 


READ  INSTRUCTIONS 
BEFORE  COMI-I.FTINO  FORM 


2.  GOVT  ACCESSION  NO.  3  RECIPIENT'S  CATAlOG  NUMBER 


im 


7  AuTHO*f«j 


8  performing  org.  report  NJMBER 


8.  CONTRACT  or  grant  n  UM8  E Rfr; 

N00014-82-K-0369 


9  PERFORMING  organization  NAME  ano  aodress 


io.  program  elEmi  •,  ■  fro„:  t»jk 

AREA  &  AO  R  K  UNiT  NUMBERS 


Mechanical  Engineering  Department  '\iRnq4-4iq 
Stevens  Institute  of  Technology  i N  9 
Hoboken,  New  Jersey  07030 _ j 


M.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 
1800  N.  Quincy  Street 
Arlington,  VA  22217 


U  MONITORING  AGENCY  NAME  &  AOORESSff/  dlltortnt  from  Controlling  OHIce)  IS.  SECURITY  CLASS,  (ol  i him  report. 

Unclassified 

See  above  _ _ _ _ _ 

IS*.  DECLASSlFlCATiON.  OOANGRAC  NO 
SCh  EDU  LE 


12.  REPORT  DATE 

15  August  1983 


13.  NUMBER  OF  PAGES 

41 


17.  DISTRIBUTION  STATEMENT  (o I  tha  abattad  antarad  In  Block  20,  II  dlllaeant  from  Raport) 


See  above 


19.  KEY  WORDS  ( Contlnua  on  ravataa  aid a  If  nacaaamry  and  Idantlly  by  block  numbar) 

Unsteady  Flow,  Compressor  Stall,  Propagating  Stall,  Rotating  Stall,] 
Laser  Doppler  Anemometer,  Laser  Doppler  Velocimeter  I 


20.  ABSTRACT  (Contlnua  on  ravataa  alda  If  nacaaamry  and  Idantlly  by  block  numbarj 

The  objectives  of  this  research  program  were  to  measure  the 
unsteady  flow  field  in  the  rotor  passage  of  an  axial  compressor 
and  map  this  velocity  field  when  operating  near  rotor  blade  stall. 
A  laser  doppler  velocimeter  system  was  developed  for  this  purpose 
and  used  in  conjunction  with  a  single  stage  compressor  designed 
expressly  to  facilitate  these  measurements.  An  extensive  data 
acquisition  system  utilizing  microcomputer  components  was  designed^ 


FORM  .... 

1  jAN  73  >4/3  EO.TiON  of  1  NOV  65  |$  OBSOLETE 
S/N  r.  102-014  -  66C  1 


UNCLASSIFIED 

security  classification  of  this  page  (iw>«n  d*i*  sn,,r. 


and  partially  assembled  and  data  processing  algorithms  were 
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